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1 .  INTRODUCTION 


There  has  been  considerable  recent  interest  in  the  use  of  nonlinear  opti¬ 
cal  materials  for  all-optical  image  and  signal  processing.*’^  The  development 
of  new  tunable-frequency  (agile)  lasers  has  also  created  demands  for  new  non¬ 
linear  materials  for  optical  harmonic  generation,  which  by  now  is  a  conven¬ 
tional  nonlinear  optical  technique.®  These  devices  or  techniques  are  ulti¬ 
mately  all  limited  by  materials.  In  this  report,  we  develop  a  theoretical 
technique  using  a  three-dimensional  anharmonic  oscillator  model  to  derive  the 
nonlinear  susceptibilities  of  crystalline  materials.  This  technique  may  have 
application  in  predicting  potential  new  nonlinear  materials.  We  derive  ex¬ 
pressions  for  the  nonlinear  susceptibility  tensor  elements  in  terms  of  the 
anharmonic  force  constants  and  derive  relationships  between  the  tensor  ele¬ 
ments  according  to  the  symmetry  of  the  32  crystal  classes. 


This  report  is  organized  as  follows.  In  section  2,  some  relevant  back¬ 
ground  material  on  nonlinear  optics  is  given.  (The  reader  is  referred  else¬ 
where**"®  for  a  more  in-depth  review.)  We  describe  our  approach  in  section 
3.  In  section  4,  the  anharmonic  oscillator  model  is  briefly  discussed.  The 
method  of  construction  of  the  anharmonic  oscillator  terms  in  the  potential 
energy  is  described  in  section  5,  with  specific  examples.  The  solution  of  the 
Lorentz  force  equation,  containing  anharmonic  terms,  is  given  in  section  6. 

r  i  \  ^  ^  311(1  '  Wo  Ho'.r, 


The  solutions  give  rise  to  x 


X'-".  We,  however,  present  the 
nonlinear  results  in  terms  of  Miller's  6.  The  final  results  for  all  the 
crystal  classes  are  collected  in  the  appendix,  where  the  invariant  polynomi¬ 
als,  resultant  potential,  and  Miller's  6  are  listed  for  each  crystal  class. 
The  results  automatically  obey  the  Kleinman  conditions.'® 


2.  BACKGROUND 

Large  optical  electric  (E)  fields  produce  a  nonlinear  response  in  the 
medium.  The  resultant  effects  are  describabie  by  the  induced  polarization 
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3)g3  +  _  _ 


(1) 


P  =  ^  +  X 


(2)ga  , 


where  x  x  ^  the  second-order  and  third-order  nonlinear  optical 

susceptibilities.  x  ^  which  we  will  refer  to  simply  as  x  is  the  linear 
susceptibility.  For  isotropic  material,  x  is  related  to  the  refractive  index 
n  by 

X  =  [n^  -  i]/4tt  (2) 

in  esu.  For  anisotropic  materials,  x  becomes  a  tensor.  We  assume  that  the 
susceptibilities  and  indices  of  refraction  are  measured  along  the  principal 
optic  axes  so  that  x^^  =  ~  1  )/4it.  We  consider  light  of  a  single  frequency 

w  and  limit  the  discussion  to  the  nonlinear  processes  of  second-harmonic  and 
third-harmonic  generation.  For  now  we  consider  second-harmonic  generation. 

The  second-order  polarization  at  2oj  is  usually  written 

=  d^^;*  ,  (3) 

1  ijk  j  k  ’ 

where  i,j,k  refer  to  the  coordinates  x,y,z  =  1,2,3  and  we  assume  summation 
over  repeated  indices.  The  are  the  second-harmonic-generation  coeffi¬ 
cients  which  are  the  components  of  a  third  rank  tensor  where  . 

Because  of  the  symmetry  of  j  and  k,  it  is  customary  to  contract  the  jk  suffix 

to  2,,  where  H  =  1  to  6  replaces  jk  =  1  1  ,  22,  33,  23  (or  32),  13  (or  31),  12 

(or  21),  respectively.  Experimentalists  generally  report  measurements  in 
terms  of  d^w  rather  than  x^^^:  however,  these  two  tensors  are  simply  related 
to  each  other  by  constants  according  to  the  particular  tensor  element; 

AV  =  At  ^  3  and  x^f  =  2d?“  for  ^  4  .  (4) 

In  terms  of  the  contracted  df^'s,  one  can  express  as 


The  symmetry  properties  of  the  crystal  classes  are  used  to  determine  the 
nonvanishing  components  and  relationships  between  elements  of  the  d^*^ 
tensor  ®  ® ^  ^  and,  similarly,  for  the  tensor,*^  as  summarized  by 

Flytzanis.’  An  immediate  consequence  of  symmetry  is  that  the  second-order 
nonlinear  coefficients  vanish  for  centres ymmetri c  media.  In  addition  to 


crystal  symmetry,  there  is  a  symmetry  relation  based  on  a  conjecture  by  Klein- 
man * ”  that,  in  a  lossless  medium,  Xijk  symmetric  under  any  permutation  of 
the  indices.  Kleinman’s  condition  for  a  general  df^  matrix  is 


16  22  23  24  14  12  , 


21  22  23  24  25  26 


31  32  33  34  35  36 


15  24  33  23  1  3  1  4 


where  we  have  only  written  the  subscripts  it. 

In  1  964,  Miller reported  that  the  quantity 


X  (2(»))x  (w)x  (w) 

ii  jj  kk 

is  nearly  a  constant  for  the  materials  that  he  investigated.  This  result, 
known  as  Miller’s  rule,  has  had  use  in  predicting  good  nonlinear  materials. 
The  rule  indicates  that  materials  with  large  linear  susceptibilities  should 
have  large  second-order  nonlinearities.  (Notice  that  the  is  defined  in 
terms  of  the  d2w  given  by  eq  (4);  otherwise  factors  of  2  appear  in  the  for 
i  >  3.) 
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To  our  knowledge,  Lax  et  al  *  ®  were  the  first  to  use  the  one-dimensional 
anharmonic  oscillator  model,  incorporated  into  a  band  model,  to  describe  the 
nonlinear  properties  of  a  solid.  Bloembergen  ‘  ®  used  the  one-dimensional 
anharmonic  oscillator  model  to  estimate  the  one-dimensional  nonlinear  force 
constant  by  considering  the  magnitude  of  the  atomic  forces  within  a  cell. 
Garrett  and  Robinson^''  extended  arguments  involving  the  one-dimensional  anhar¬ 
monic  oscillator  to  make  order-of-magnitude  calculations  of  Miller's  6. 
Yariv‘®  provides  a  concise  review  of  this  work.  Adler, in  his  discussion  of 
nonlinear  optical  frequency  polarization  in  a  dielectric,  introduces  constants 
in  the  force  equations  of  a  three-dimensional  harmonic  oscillator  model; 
however,  he  does  not  investigate  the  properties  of  the  different 

crystal  classes.  Robinson^”  discusses  the  three-dimensional  harmonic  oscilla¬ 
tor  model  but  does  not  investigate  the  consequence  of  making  the  potential 
energy  invariant  for  the  particular  crystal  classes.  He  does,  however,  point 
out  the  relationship  of  the  constants  in  the  potential  if  the  potential  satis¬ 
fies  Laplace's  equation.  We  have  not  found  any  detailed  application  of  the 
anharmonic  oscillator  model  in  three  dimensions  to  nonlinear  optics. 

3.  APPROACH 

In  this  report  we  introduce  anharmonic  terms  into  the  potential  energy  of 
a  three-dimensional  oscillator  model  such  that  each  term  is  invariant  under 
the  operations  of  the  point  group  of  a  particular  crystal  class. The  32 
point  groups  are  listed  in  table  1.  Invariant  polynomials  of  order  two  are 
sufficient  for  the  linear  susceptibility,  x:  third-degree  polynomials  are 
required  for  x^^^  and  fourth-degree  polynomials  are  used  for  f’on  the 

lowest  symmetry  crystal  class  (triclinic,  (1))  the  number  of  independent 
polynomials,  'Nj^,  of  degree  n,  is  given  by  the  number  of  terms  in  the  expansion 
of  (x  +  y  +  z)n.  For  a  given  n,  there  are  =  (n  +  1)(n  +  2)/2  polynomials. 
For  the  lower  crystal  classes  (triclinic,  monoclinic,  and  orthorhombic — 1 
through  8),^  we  assume,  with  no  loss  of  generality  (but  a  great  deal  less 
algebra),  that  the  coordinate  system  x,y,z  is  parallel  to  the  principal  optic 
axes  of  the  crystal.  Using  the  potential  energy  in  the  Lorentz  force  equa¬ 


ls 
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tThe  ordering  of  the  crystal  classes  is  as  given  by  Koster  et  al  (ref  21). 


TABLE  1  .  CLASSIFICATION  OF  THE  32  POINT  GROUPS 


System 

Unit  cell 

Point 

group 

number 

Symmetry 

Number  of 

symmetry 

elements 

TrLclinio 

a  * 

b 

»  c 

1 

C,  (1) 

1 

a  ■ 

S 

»  y 

2 

(T) 

2 

Monoclinic 

a  * 

b 

*  0 

3 

Cj  (2) 

2 

a  ~ 

Y 

=  tt/2  »  6 

1) 

Cs 

2 

5 

C2^  (2/m) 

9 

Orthorhombic 

a  * 

b 

*  c 

6 

Dj  (222) 

9 

a  = 

6 

=  Y  =  ti/2 

7 

Cjy  (mm2) 

9 

8 

Djf,  (mmra) 

8 

Tetragonal 

a  = 

b 

»  C 

9 

C4  (9) 

9 

a  = 

B 

=  Y  =  Tt/2 

to 

S4  (9) 

9 

1 1 

Cuf,  (Vm) 

8 

12 

(  922) 

8 

13 

(9mm) 

8 

19 

Djj  (92m) 

8 

15 

(9/mmm) 

16 

Rhombohedral 

a  = 

b 

*  C 

16 

C3  (3) 

3 

(trigonal ) 

a  = 

6 

=  Y  <  2Tr/3  •  ii/2 

17 

Csi.Sfe  (3) 

6 

18 

D3  (32) 

6 

19 

C3V  (3ni) 

6 

20 

D3J  (3m) 

12 

Hexagonal 

a  * 

b 

»  c 

21 

Cg  (6) 

6 

a  » 

B 

=  ii/2,  Y  =  2it/3 

22 

C3h  (6) 

6 

23 

Cbh 

1  2 

29 

Dj  (622) 

12 

25 

Cgy  (6mra) 

12 

26 

03^  (6m2) 

12 

) 

27 

*^60  (^^mmm) 

29 

Cubic 

a  * 

b 

»  c 

28 

T  (23) 

12 

a  = 

B 

=  Y  =  11/2 

29 

(ra3) 

29 

30 

0  (932) 

29 

31 

T(j  (93m) 

29 

32 

0^  (m3m) 

98 

tion,  the  susceptibilities  x»  X^  ^  derived  for  the  32  crystal 
classes.  We  show  that  all  the  susceptibilities  thus  derived  obey  the  Klein- 
man's  conditions'®  given  in  equation  (6).  In  fact,  for  any  crystal  class,  the 
number  of  independent  third-  or  fourth-degree  polynomials  is  equal  to  the 
number  of  independent  df^  or  d?^,  so  that  each  corresponds  to  a  particu¬ 
lar  constant  in  the  enharmonic  potential.  No  attempt  is  made  here  to  derive 
the  constants  in  the  potential  energy  from  more  fundamental  theory.  That  is, 
each  constant  in  the  potential  energy  allowed  by  symmetry  is  assumed  phenome¬ 
nological  and,  at  this  point,  is  determined  by  fitting  experimental  data.  In 
particular,  the  constants  in  the  quadratic  terms  of  the  potential  energy,  the 
harmonic  oscillator,  can  be  determined  from  the  measured  index  of  refraction 
at  various  wavelengths  by  using  the  Sellmeyer  or  Sellmeier  equations . ^ ^ 


The  6ijk  of  Miller's  rule  are  given  simply  in  terms  of  the  constants  (6^) 
of  the  potential  energy  with  n  =  3  so  that  these  constants  can  be  related 

to  experimental  data.  The  extension  of  Miller's  rule  to  is  straightfor¬ 
ward  for  those  crystal  classes  where  =  0  (6i  =  0),  in  which  case  the 

simply  related  to  the  constants  (Y^)  of  the  fourth-degree  terms  in 
potential  energy.  The  extension  of  Miller's  rule  for  to  those  crystal 

classes  where  6?*^^  *  0,  which  is  generally  not  simple,  is  also  given. 


4.  ANHARMONIC  OSCILLATOR  MODEL 

We  assume  that  the  nonlinear  polarization  is  electronic  in  origin  and 
start  with  the  Lorentz  model  which  has  been  used  to  describe  the  linear  re¬ 
sponse  of  the  electrons  in  a  solid  to  an  electric  field. In  this  model,  the 
equation  of  motion  for  an  electron  is 

mp  =  -eE  -  VD  ,  (8) 

where  m  is  the  mass  of  the  electron,  e  is  the  electronic  charge,  r  is  the 
position  of  the  electron,  E  is  the  applied  electric  field  seen  by  the  elec¬ 
tron,  and  U  is  the  potential  of  the  electron  in  the  solid.  We  will  ignore 
losses  which  could  be  included  by  a  damping  term  -mfr  on  the  right-hand  side 
of  equation  (8)  where  f  is  the  damping  constant.  In  an  isotropic  solid,  the 
potential  for  a  harmonic  oscillator  model  is  chosen  as 

U  =  ^  mwQ(x^  +  y2  +  z^)  ,  (9) 


where  Wq  is  a  representative  resonance  frequency  of  the  solid, 
tial  given  in  equation  (9)  is  used  in  equation  (8),  we  obtain 

mr  =  -eE  -  mwgr  , 


If  the  poten- 


^D.  A.  Kleinman,  Nonlinear  Dielectric  Polarization  in  Optical  Media,  Phys.  Rev.  126  (1962), 


^^M.  Born  and  E.  Wolf,  Principles  of  Optics,  Pergamon  Press,  New  York  (1964),  p  97. 

F.  J.  Bottcher  and  P.  Bordewijk,  Theory  of  Electric  Polarization,  Vol  II,  Elsevier,  New 
York  (1978)  ,  p  288. 

0.  Jackson,  Classical  Electrodynamics,  2nd  Edition,  McGraw-Hill,  New  York  (1975),  p  285. 


and  if  a  time  dependence  of  cos  wt  is  chosen  for  E,  we  obtain 


2~ri/N*  (11) 

Wq  -  (Jj  Dilo) 

where  D((jj)  =  -  uj^j/e.  The  polarization,  P,  is  related  to  the  displace¬ 

ment,  r,  by 

P  =  -Ner  ,  (12) 

where  N  is  the  number  of  electrons  per  unit  volume  that  contribute  to  P. 
Substituting  (11)  into  (12)  and  using  P  =  xE,  we  obtain 


X(w) 


Ne 

b(w) 


(13) 


For  large  displacement  from  the  electronic  equilibrium  position,  the 
anharmonicity  of  the  electron  oscillators  must  be  taken  into  account  as  done 
by  Bloembergen ' ®  for  the  one-dimensional  anharmonic  oscillator  model.  In  the 
lowest  order  nonlinear  approximation,  the  anharmonicity  is  introduced  through 
a  potential  energy  term  proportional  to  the  cube  of  the  displacement — 
[“(1/3)Sx^]  —  where  we  use  6  for  the  nonlinear  force  constant.  This  term  gives 
rise  to  second-harmonic  generation  through  the  nonlinear  polarization  term 


=  X^^^E5j(w)E^(6o)  , 

where  the  second-order  nonlinear  susceptibility  can  be  shown  to  be 


(1U) 


(15) 


(equation  (1-14)  of  Bloembergen * ® ) . 


5.  ANHARMONIC  POTENTIAL 

5 . 1  Invariant  Polynomials 

By  far  the  most  difficult  part  of  the  work  reported  here  is  the 
construction  of  the  anharmonic  terms  in  the  potential  energy.  Unfortunately, 
we  have  not  found  a  universal  method  of  constructing  these  polynomials  but 
have  resorted  to  a  number  of  devices  to  achieve  our  goal.  The  most  useful 
technique  is  what  might  be  called  the  replacement  method.  That  is,  for  a 
given  operator,  0,  of  a  group  we  have 
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N. 


V. 


Bloembergen , 
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0(x,y,z)  ■*  (some  permutation  of  x,y,z)  (16) 

for  each  operation  0  of  the  group.  For  example,  the  operator  C|^,  which  repre¬ 
sents  a  90°  rotation  about  z,  can  be  written 

Cn(x,y ,z)  -*■  (y ,-x,z)  ,  (17) 

The  32  point  groups  are  listed  in  table  1.  For  the  operations,  we 
follow  the  conventions  of  Koster  et  al,^’  where  an  explicit  description  of 

each  of  the  operations  of  the  point  groups  is  given.  Each  polynomial  and 

combination  of  polynomials  can  be  tested  for  invariance,  and  proper  combina¬ 
tions  can  be  selected.  This  method  works  quite  well  for  all  the  groups  except 
for  those  with  a  threefold  or  sixfold  rotation.  For  other  conventions,  such 
as  IRE,  the  reader  can  extend  the  methods  used  here  in  a  straightforward 
manner . 

It  is  convenient  to  list  explicitly  the  polynomials  of  various  degree 
in  tabular  form.  The  polynomials  of  order  2,  3,  and  are  given  in  table  2, 

which  also  identifies  the  labels  of  the  coefficients  of  the  polynomials  of 

order  3  and  4,  g  and  Y,  respectively. 

5. 2  Examples 

5.1.2  Crystal  Class  10, 

Crystal  class  is  a  cyclic  group  containing  the  single  generator 
Si|  which  gives 


S4(x,y,z)  =  (y,-x,-z)  , 

S^(x,y,z)  =  (-x,-y,+z)  ,  (18) 

S’(x,y,z)  =  (-y,x,-z)  . 

We  see  immediately  that  is  invariant  and  x^  +  is  also  Invariant. 

(2) 

Using  table  2a,  we  can  now  write  the  term  in  the  potential,  as 

+  y^]  +  4  ™ 

2  x  X  2  z  z 

The  m^^  and  m^  are  introduced  as  "effective  masses"  which  can  be  determined  by 
fitting  the  linear  susceptibilities  to  the  measured  indices  of 
refraction .22,23 


21 
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Groups,  MIT,  Cambridge,  MA  (1963). 
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C.  F.  J.  Bottcher  and  P.  Bordewijk,  Theory  of  Electric  Polarization,  Vol  II,  Elsevier,  New 
York  (  1978) ,  p  288. 


12 


(a)  n  =  2  where  *  T  a.Q. 

“11 


r-n 


«i 

(2) 

i 

2 

m  u) 

X  X 

1  2 

2 

2 

m  0) 

y  y 

1  2 

2  y 

2 

m  cij 

z  z 

1  2 

2  " 

n  =  3  where  = 

6. 

1 

q(3) 

1 

®2 

- 

3  ^ 

®3 

2 

X  y 

«5 

2 

X  z 

®6 

2 

y  X 

^7 

) 

2 

y  z 

®8 

2 

Z  X 

®9 

2 

z  y 

^0 

2xyz 

1  1 


u(n) 

n  =  ^ 

where 

1 

"1 

1  4 

4 

"2 

1  4 

4  y 

^3 

1  4 

4  ^ 

x^y 

^5 

x^z 

y^x 

y'z 

^8 

z3x 

^9 

z'y 

^0 

|xV 

^11 

322 

2  "  " 

"12 

322 

2  y  " 

^13 

3x  yz 

^14 

,  2 

3y  xz 

^15 

,  2 

3z  xy 

'V 


tv-jO<k.W.-C' V-.' ' 


We  find  that  xyz  and  z(x^  -  y^)  are  the  only  invariant  polynomials 
of  order  three.  Using  table  2b  we  have  Bj  =  -65  and  consequently, 

=  p^(x^  -  y^)z  +  26Qxyz  (20) 

and  all  other  3^  =  0. 

For  the  fourth-degree  terms  in  the  potential  energy,  we  see  immedi¬ 
ately  that  z‘*  is  invariant.  Since  x^  +  y^  is  invariant,  (x^  +  y^)^  is  also; 
however,  x”  +  y"  is  invariant,  so  we  additionally  have  x^y^  and  z^(x^  +  yM. 
The  last  remaining  fourth-degree  polynomial  is  xy(x^  -  y^).  Thus,  the  invar¬ 
iant  polynomials  of  order  are 

x-  +  y",  z\  xy(x^  -  y^),  x^y^,  z^(x^  +  y^)  .  (21) 

From  table  2c  we  have 

^2  ~  '*'1  >  ~  '^12  ~ 

(22) 

Y5  =  Yy  =  Yg  =  Yg  =  Y^3  =  Y^^  =  Y,5  =  0 

with  the  resulting  potential  energy 

1  y-)  ^  ^  ygZ-  Y^xy[x"  -  y") 

+  —  Y  x^V^  +  —  Y  z*(x^  +  V^)  . 

2  10  ^  2  11  ^  y  i  - 

The  results  given  in  equations  (19),  (20),  and  (23)  express  the  potential 

energy  expanded  through  polynomials  of  fourth  order. 

5.2.2  Crystal  Class  16,  Cg 

The  crystal  class  Cg  is  a  cyclic  group  containing  the  single  gener¬ 
ator  Cg  which  gives 

Cg(x,y,z)  -►  (x  cos  4)  +  y  sin  (p,  -x  sin  b  +  y  cos  (p,  z) 

(2H) 

C^(x,y,z)  -►  (x  cos  2<J)  +  y  sin  2(p,  -x  sin  2(p  +  y  cos  241,  z) 

with  4>  =  211/3.  Since  z  is  invariant,  any  power  of  z  is  invariant.  Note  that 
x^  +  y^  is  invariant,  so  that  the  potential  of  second  order  can  be  written 
(table  2a) 

=  -5  '''  u)^(x^  +  y^)  +  4  •  (25) 

2  x  X  2  z  z 

For  the  third-degree  polynomials,  we  have  immediately  z’  and  z(x^  +  y^).  To 
obtain  the  other  polynomials  using  the  replacement  method,  the  amount  of 


algebra  is  prohibitive,  and  we  resort  to  methods  used  in  the  theory  of  para¬ 
magnetic  ions  in  crystal ^  From  these  references  and  using  Koster  et 
al^‘  (page  53).  we  find  that  the  spherical  harmonics  Y^g  and  Y^  ^2  ^re  invar¬ 
iant  under  the  operation  C^.  Tables  of  the  spherical  harmonics'in  Cartesian 
form  are  given  by  Ballhausen*®  and  Prather. Leaving  out  the  normalization 
constants,  they  are 


Y30  '  2z^  -  3z(x^  +  y^)  , 

Real  Y22  ~  -  3xy^  ,  (26) 

Imaginary  Y^^  ~  3x^y  -  y^  . 

We  already  have  found  the  two  polynomials  z’  and  z(x^  +  y^)  of  Y^g,  but  the 

two  contained  in  Y^^  are  new.  No  other  combinations  can  be  found.  We  now 

have  all  the  polynomials  of  third  order:  z\  x^  -  3xy^,  y’  -  3x^y,  and 

z(x*  +  y^).  From  table  2b  we  have  =  -B-]  ,  B4  =  “62’  ^7  "  ^5’  ^  ^9  " 

8g  =  0.  Thus,  we  have 

^  6^(x'  -  3xy^)  82(7"  -  3x^y]  +  ^  B^z®  +  B^ztx^  +  y^]  (27) 

for  the  potential  energy  of  third  order. 

The  terms  of  fourth  order  in  the  potential  are  obtained  using  the 
previous  methods  and  are 

y^]^  ,  z-  ,  z(x’  -  3xy^]  ,  2(y*  -  3x^y)  ,  and  z^[x^  +  y*)  .  (28) 

These  results  can  be  checked  by  examining  the  expressions  for  Y^g  and  Y^^ 
(real  and  imaginary  parts)  which  if  done  shows  that  there  are  no  more 

polynomials  of  fourth  order.  The  constants  in  the  potential  energy  involving 
the  fourth-order  polynomials  given  in  equation  (28)  can  be  determined  using 
table  2c: 

^2  "  ’  ^10  "  3  "  ”^5  ’  ^13  "  ~^7  ’  "^12  "  ^11 

and  (29) 

Y4  =  ^6  =  ^8  =  ''9  =  ''15  =  0  ' 

The  fourth-ord-^r  potential  energy  is  then  given  by 


21 
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(30) 


\  ^  Y^z(x^  -  3xy^] 

+  Y^z(y’  -  3x"y)  +  |  Y^^z^(x^  +  y^]  . 


The  invariant  polynomials  for  all  the  32  crystal  classes  can  be  constructed 
using  the  above  techniques.  In  appendix  A,  the  invariant  polynomials  and 
potential  energy  expanded  through  order  4  are  explicitly  written  for  each 

^  in 


crystal  class.  The  results  are  given  in  summary  form  for  U 
for  in  table  4. 


table  3  and 


6.  EQUATION  OF  MOTION  AND  SOLUTION 


To  simplify  the  solution  of  the  equation  of  motion,  equation  (8),  we 
assume  that  the  material  is  lossless  in  the  region  of  interest.  This  reduces 
the  algebra  considerably.  We  use  a  perturbation  technique  which  is  essen¬ 
tially  an  expansion  of  the  displacements  in  powers  of  the  electric  field.  To 
illustrate  the  technique  we  shall  set  up  the  equations  and  their  solution  for 
the  C-j  crystal  class.  Since  the  crystal  class  has  no  symmetry,  the  results 
for  any  crystal  class  of  higher  symmetry  can  be  obtained  by  using  the  methods 
of  the  previous  section  along  with  the  results  of  tables  3  and  4. 


From  appendix  A  the  potential  energy  for  the  crystal  class  C-j  is  given  as 


(31 ) 


+  B^x^z  +  S^xy^  +  B^y^z  +  BgXZ^  +  B^yz*  +  2BQXyz  +  ^  Y^x" 

+  ^  Y^y'*  +  ’Ij  Y^z**  +  Y^x*y  +  YgX^z  +  YgXy®  +  Y^y’z  +  YgXz’ 

*  V^’  *  I  to’'”''  *  I  *  I  *  3'f,3x‘yz 

+  3Y^;^xy^z  +  3Y^^xyz^  . 

The  components  of  the  equation  of  motion  are  obtained  by  using  the  above 
expression  in  equation  (8)  to  give 


mvX  = 


eEx  -  m^ioxx  -  Bix^  -  2B4xy  -  2B5XZ  -  ~  632  '  2Boy^ 


-  Y^x-^  -  3Y,x^y  -  3Y5X^z  -  Y^y^  -  YgZ^  -  3Y^oXy^ 

-  3Y^^xz2  _  bY^^xyz  -  3Y^^jy^z  -  3Y^^yz^ 


(32) 


with  similar  equations  for  y  and  z. 


We  assume  that  the  field  and  that  x  =  x^6  +  X26^  +  XgB^,  with  a 
similar  expression  for  the  other  two  components,  and  where  6  is  such  that  at 
the  end  we  set  6^1. 
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TABLE  3.  FOR  ALL  CRYSTAL  CLASSES'^ 


No.  Symmetry  6 


“2h 


”2v 


‘"2h 


'Mh 


'4v 


^2d 


'3v 

“3d 


'3h 


“6h 


■"Sv 


6h 


Oh 


1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-6, 


0 

0 

0 

-82 

0 

0 

0 

-82 

0 

0 

0 

0 

0 

0 


1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-8, 

0 

0 

-8, 

0 

0 

-8, 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


1  1 
0  0 
0  1 
1  0 
0  0 
0  1 
0  0 
0  0 
0  0 
0  1 
0  0 
0  0 
0  0 
0  1 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  0 
0  1 
0  0 
0  0 
0  1 
0  0 


i"  or  "0"  in  the  tebie  represents  the  multiplier  of  B  •  Other  entries 
"-Bj";  represent  the  replacement  term  for 
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TABLE  H.  FOR  ALL  CRYSTAL  CLASSES'* 


No.  Symmetry 


'11  '12  '13  14 


11111 


>,/3  1 

0 

v,/3  1 

"11 

0 

Y,/3  1 

'^11 

0 

1 

-< 

0 

Y,/3  1 

''n 

0 

0 

Y,/3  1 

0 

-''5 

0 

Y,/3  1 

^1 

0 

0 

0 

Y,/3  1 

'n 

0 

0 

0 

Y,/3  1 

0 

0 

0 

Y/3  1 

''n 

0 

0 

0 

Y,/3  1 

''n 

0 

0 

0 

Y,/3  1 

0 

0 

0 

Y,/3  1 

'"n 

0 

0 

0 

0  0 
0  0 


0 


0 


0 


Using  these 

expressions  in  equation 

(32)  and  equating  powers  of 

6 ,  we 

obtain 

mxXi  =  -eEx  - 

'"x‘^x^1  > 

myyi  =  -eEy  - 

my'jjyyi  , 

(33) 

^Z^l  ~  '^^2 

^z'^z^l  > 

for  the  terms  of  first  order  in  6, 

V2  =  -\<^2  '  ■  265^1^1  ~  Vi'  ■  ®8^l'  ■  ^Vl^  ^34) 

for  the  terms  of  second  order  in  6,  and 


™x^3  =  ■  261X1X2  -  26ii(xiy2  ^  X2yi  j  -  2B^[x^Z2  *  X2Z1  ) 

-  3Y5X^^z^  -  Y^y^^  -  Ygz’  -  3Y^QX^y^  -  3Y^^x^z^^ 


(35) 


-  6Y^3X^y,z^  -  3Y^^y^^z^  -  SY^^y^z,^ 

for  the  terms  of  third  order  in  6  where,  for  simplicity,  we  have  written  only 
the  x-components . 


The  results  given  in  equations  (33).,  (3^),  and  (35)  are  sufficient  to 
obtain  the  susceptibilities  x>  X  >  ^nd  x  We  make  use  of  equation  (12) 

generalized  to 

=  -Nex  (36) 

X  n 

with  n  =  1,  2,  3  and  similar  expressions  for  y  and  z. 

We  start  by  solving  the  first-order  equations  given  in  equation  (33). 
Since  we  are  Ignoring  losses,  we  can  assume  that  any  time-dependent  quantity 
varies  as  co3((i)t  -  <()) ,  so  that  the  solutions  are 

^1  ~  ~^X^^X  ’ 

y,  =  "V'^y  ’  ^37) 

^1  "  "^z^^z  » 

and  consequently,  we  obtain  the  components  of  the  linear  susceptibility  tensor 

Xii  =  Ne/D^  ,  (38) 

where  =  D^(u))  =  m.  (u?  -  oi^j/e.  We  will  give  the  argument  of  explicitly 
only  when  it  is  2u  or  3u. 
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Using  the  result  of  equation  (37)  in  equation  (3^),  we  have 

6nE^  26 „E  E 


2  eD  (2cj)D^  eD  (2a))D^  eD  (2(a))D^  eD  (2oj)D  D 
X  X  X  y  X  2  x  xy 

26cE  E  26„E  E 

5  X  z  0  y  z 

eD~r2(J)D  D  ~  eD  (2w)D  D 
X  X  z  X  y  z 


We  can  rewrite  equation  (5)  as 


I  d->c 


where  \f ^  are  given  in  terms  of  the  electric-field  components  according  to 
equation  (5).  Using  equations  (36),  (38),  and  (^0),  we  obtain 


^2  "  ^^6’^xx^^‘"^>^yy)<yy  * 


where  g  =  1/N^e’.  Comparing  the  results  given  in  equation  (41)  with  equation 
( 7) ,  we  have 

6ll  =  ,  6^2  "  “  ^8®  ’  “^IM  "  ^0^  ’  '^15  °  ’ 


. 

Similarly,  the  solutions  for  y2  and  Z2  give  rise  to  62^  and  63^*  We  can 
write  the  resultant  (Sf^  in  matrix  form  as 


.20) _ 1 

,,2  3 
N  e 


B5  S4 


62  Sg  8^  6q  Bg  . 


®7  ^3  ®8  ®0 


Note  that  satisfies  the  Kleinman  conditions  given  in  equation  (6) 


The  solutions  for  the  third-order  displacement  obtained  by  using 

the  results  of  (37)  and  (39)  in  (35).  It  can  be  seen  by  inspection  of  equa¬ 
tion  (35)  that  we  can  express  the  solution  as  =  x^  +  x'^  where  x^  ^  0  as 
-♦  0  and  xU  ->0  as  -►  0.  For  simplicity  of  demonstration,  we  will  write 
only  the  portion  of  x^  containing  the  Y  terms,  which  is 


. .  V)  V)  3Y,E»E 

eD  (3a.)x5  . 

X  y  z  X  y 


3Y,  ,,E^E  3Y,,E^E  3Y^E^E 

,14yz,11zx.5xz 

D^D  ^  D  ^  D^D 

Y  Z  X  Z  X  z 


3Y, „E  E^  3Y, ^E  E^  6Y, ^E  E  E 
10  X  y  ^  15  y  z  ^  13  X  y  z 

D  D^D  D  D  D 

X  y  y  z  X  y  z 


The  solution  for  x,  gives  rise  to  P' 
3  °  > 


The  third-order  polarization  can  be  expressed  as  Pj  =  Xi j ;  uuw- 
ever,  in  analogy  to  the  case  for  P^***,  we  use  the  reduced  oasis  for  P^'*'. 
Following  the  convention  of  Maker  and  Terhune,^®  we  write  P^“  =  Xim  ^m  where 

subscripts  jkH  are  rep'aced  by  m  and  EjE|^Ej^  is  replaced  by  Uj^^,  as  shown  in 

table  5.  Again  in  analogy  with  second-harmonic  generation,  we  can  express 

p?^  =  ^3*^0  (441 

where  d3j^  are  the  third-harmonic  generation  coefficients  and 

Qm  =  Ujn  for  1  1  m  <  3  , 

Qjj,  =  3U^  for  >3  <  m  <  9  ,  (45) 

Qg  =  6Uq  for  m  =  0  . 

TABLE  5.  RELATIONSHIP  BETWEEN  P?“  =  x-?»E.E,  E„  and 

1  '^ijkl  j  k  a  *im  m 


,E.E|.E|,;  hgw- 


P.  0.  Maker  and  R.  W.  Terhune.  Study  of  Optical  Effects  due  to  an  Induced  Polarization  Third 
Order  in  the  Elect-1c  Field  Strength,  Phys.  Rev.  137A  (1965),  801. 


Using  the  results  given  in  table  5  along  with  equations  (4i4)  and  (45),  we 
can  write  the  Kleinman  conditions  for  as 
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Generalization  of  Miller's  rule  to  the  case  of  the  third -harmonic  generation 
gives^®’  ° 


ijkJ,  ijk2,  ‘■''11  ''jj'' 

Using  equations  (42),  (43),  and  (44),  we  obtain 


d3»  . 

ofj  ■ 

d3w  - 

di3 


Xxx^ Xxx  ’ 

Xyy  f 

-bY4Xxx(3a))x^,Xyy  .  etc. 


where  b  =  l/N^e**. 

Note  that  an  advantage  of  using  the  d^^  rather  than  x^^^  ie  that  the 
factors  of  3  and  6  which  occur  in  equation  (43)  are  eliminated. 

Comparing  the  terms  in  (48)  with  equation  (4?)  gives 
5 1 1  =  "t)Yi  ,  1 2  ~  ’  *^13  ~  ’  ^^14  ”  >  *^15  ~  ~^''^14  ’ 

(49) 

'^16  "  ”^''^11  ’  '^17  "  '^’'^5  ’  '^18  “  "'^^10  ’  '^19  ^  ”^^15  ’  "^10  "  “*^^13  ■ 

When  the  results  for  the  y''  and  z''  solutions  are  included,  one  can  express 
63^  in  matrix  form 


J.  J.  Wynne  and  G.  D.  Boyd,  Study  of  Optical  Difference  Mixing  in  Ge  and  Si  using  a  CO2  Gas 
Laser,  Appl.  Phys.  Lett.  12  (1968),  191. 

C.  Wang,  Empirical  Relation  Between  the  Linear  and  the  Third-Order  Nonlinear  Optical 
Susceptibilities,  Phys.  Rev  B2  (1970),  2045. 
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The  reader  should  recall  that  this  solution  ignores  the  B  contribution.  The 
complete  solution  is  given  in  the  appendix  and  obeys  the  Kleinman  condition 
given  in  equation  (M6). 


7.  CONCLUSIONS 

We  have  derived  the  nonvanishing  coefficients  for  optical  doubling  and 
tripling.  The  derivation  is  based  on  the  Lorentz  model  where  the  linear  re¬ 
sponse  of  the  bound  electrons  to  an  incident  electric  field  is  described  by  a 
classical  harmonic  oscillator.  Anharmonic  terms  (i.e.,  terms  of  higher  order 
than  quadratic  in  the  potential)  are  added  to  describe  the  nonlinear  response 
of  the  medium.  These  anharmonic  terms  are  subject  to  the  restriction  that  the 
potential  remain  invariant  under  the  operations  of  each  particular  crystal 
class.  The  nonlinear  optical  constants  dfjj^  and  were  obtained,  and  the 

results  were  expressed  in  terms  of  Miller’s  6's.  Explicit  expressions  for  the 
6's  are  given  in  terms  of  the  coefficients  of  the  third-  and  fourth-degree 
polynomials  of  the  potential.  These  coefficients  were  considered  phenomeno¬ 
logical,  and  no  attempt  was  made  to  derive  their  values  from  a  theoretical 
model . 

Additional  assumptions  can  be  used  to  place  further  restrictions  on  the 
coefficients  of  the  anharmonic  terms  in  the  potential.  If  we  assume  that  the 
electrons  contributing  to  the  nonlinear  susceptibilities  are  in  a  region 
relatively  free  of  charge,  then  the  potential  obeys  Laplace's  equation,  that 
is,  =  0  for  each  n.  This  condition  results  in  the  following  relations 

between  the  coefficients  of  the  potential  for  n  =  3: 

6^  +  65  +  8g  =  0  ,  62  Bij  +  Sg  =  0  ,  and  Bg  +  Bg  +  By  =  0 

For  the  crystal  classes  Cg  and  C6v.  the  condition 

Bg  Bg  By  “  0  , 

when  used  in  our  results  for  Miller's  6fk,  requires  that  263“  +  633  =  0.  This 
result  agrees  with  experimental  data  (within  reported  error)  tabulated  by 
Singh®  for  five  of  the  seven  values  for  the  crystal  class  Similarly,  for 

n  =  4,  additional  restrictions  can  be  found  for  the  Y^^'s. 

^S.  Singh,  Non-linear  Optical  Materials,  in  Handbook  of  Lasers,  ed.  by  R.  J.  Pressley,  The 
Chemical  Rubber  Co.  (1971). 


A  number  of  authors  have  used  quantum  mechanics  to  derive  expressions  for 
the  nonlinear  coefficients  (see  Flytzanis,^  for  example);  however,  the  useful¬ 
ness  of  this  approach  has  been  limited  by  the  inability  to  evaluate  the  re¬ 
quired  matrix  elements.  Although  the  classical  anharmonic  oscillator  model 
has  been  applied  to  nonlinear  optics  by  a  number  of  authors,  to  our  knowledge 
none  has  derived  explicit  expressions  for  the  nonlinear  coefficients  in  terms 
of  the  anharmonic  coefficients  for  each  crystal  class.  It  should  be  possible 
to  explicitly  calculate  the  magnitude  of  and  ,  which  would  be  useful  in 
predicting  new  nonlinear  materials.  Work  is  in  progress*  in  calculation  of 
these  coefficients  using  a  point-charge  model  of  crystal  fields  in  an  ionic 
solid. 
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APPENDIX  A. --SUMMARY  OF  RESULTS  FOR  ALL  CRYSTAL  CLASSES 


The  invariant  polynomials,  resultant  potential,  linear  susceptibilities, 
and  Miller's  and  are  given  according  to  crystal  class.  From  the 
second-harmonic-generation  coefficients  d^*'*  can  be  constructed  using  equation 
(7)  which  we  can  rewrite  as 


Note  that  the  appendix  makes  use  of  the  contracted  indices  where  the  jk  suffix 
of  is  replaced  by  Si,  where  £,  =  1,  2,  3.  5,  6  replaces  jk  =  11,  22,  33i 
23  (or  32),  13  (or  31),  12  (or  21),  respectively,  x  ^  can  be  determined  from 
d^“  according  to  equation  (4)  in  the  body  of  the  report,  repeated  here: 


(2) 

'iil 


.2uj 


for  Si  <  3  and 


(2) 

’^i£ 


■  "“u 


for  i  >  4 


Similarly,  the  third-harmonic-generation  coefficient  d^“  can  be  obtained 
from  the  according  to  equation  (47),  which  is  rewritten  as 
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3u) 

ijkH 


Again  in  the  appendix  we  use  the  contracted  indices  where  m  replaces  the 
jkil  indices  according  to  table  5  in  the  text.  x  ^  can  then  be  determined 
fr  m  d^'^  according  to  x[m  "  ‘^im  ^  3,  3d^^  for  3  <  m  $  9,  and  6dj^^  for 
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Invariant  polynomials  of  fourth  order: 
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